Excitations in correlated superfluids near a continuous transition into a supersolid 
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We study a superfluid on a lattice close to a transition into a supersolid phase and show that a 
uniform superflow in the homogeneous superfluid can drive the roton gap to zero. This leads to 
supersolid order around the vortex core in the superfluid, with the size of the modulated pattern 
around the core being related to the bulk superfluid density and roton gap. We also study the 
electronic tunneling density of states for a uniform superconductor near a phase transition into a 
supersolid phase. Implications are considered for strongly correlated superconductors. 

PACS numbers: 71.27.+a, 05.30.Jp, 03.75.Lm 
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The theoretical idea that 4 He could become a super- 
solid at low temperature has motivated experimental 
studies of defect excitations in the solid phase [2j , as well 
as a theoretical investigation of vortices in a superfluid 
near a first-order transition to a supersolid 0], as indi- 
rect windows into the proposed supersolid phase. Recent 
experimental hints [4] for a non-classical rotational iner- 
tia in solid 4 He have revived interest in this field. The 
copper oxide superconductors (SC) present a different 
situation, where the primary interest is in the uniform 
superconducting state, but it has long been recognized 
that the underdoped regime is plagued by various com- 
peting phases || involving spin, charge, and current or- 
dering. Such competing phases can significantly affect 
the excitations in the superconducting state if there is a 
continuous transition between the uniform SC and these 
ordered phases. This is viewed as one possible cause 
for the anomalous excitation spectra of the underdoped 
cuprates. However, continuous phase transitions between 
phases with different order parameters are unusual and 
have only recently begun to be explored 0. Further, 
recent tunneling experiments Q in highly underdoped 
Caa-zNa^CuC^Ck, (T c « 10- 20K) suggest that the 
ground state in this regime may be a supersolid rather 
than an insulator. Motivated by this, we consider here 
a continuous transition between a superfluid (supercon- 
ductor) and a supersolid, on a lattice, and explore some 
qualitative consequences for excitations in the uniform 
superfluid (superconducting) state near such a transition. 

The following arc our main results, (i) A continuous 
transition into a commensurate supersolid phase is known 
to arise, with increasing interactions, from condensation 
of rotons when the roton gap vanishes. Here, we show 
that alternatively a uniform current flow in the superfluid 
can also induce supersolid order by driving the roton gap 
to zero, (ii) This current-driven collapse of the roton 
gap results in a supersolid pattern emerging around the 
vortex core in the uniform superfluid as shown in Fig. ^ 
The length scale over which this modulation is significant 
is R v ~ D s /E Ia t, where D s is the bulk superfluid stiff- 
ness and E TOt is the roton gap in the homogeneous super- 
fluid, (iii) We use microscopic calculations and Ginzburg- 




FIG. 1: (color online) Profiles of the modulated (p, IV'qI) an d 
uniform (m, \ipo\) components of the density/superfluid order 
parameters around a vortex in a superfluid near a transition 
to a supersolid. Inset: Density modulations near the core. 

Landau (GL) theory to obtain the long distance profiles 
of the modulated and uniform components of the den- 
sity and superfluid order parameter around the vortex, 
as shown in Fig. ^ (iv) Finally, we ask whether tun- 
neling into a uniform SC near a transition into a super- 
solid phase can probe the low energy roton excitations. 
Within a slave boson approach, we show that a tunnel- 
ing electron can excite the bosonic condensate modes in 
the SC. In SCs with a small superfluid density, this leads 
to inelastic secondary peaks in the tunneling spectrum 
as shown in Fig. [S] Coherence factors however suppress 
contributions from wavevectors corresponding to the ro- 
ton minimum. We conclude with possible implications 
for correlated SCs such as the doped cuprate materials. 

Microscopic model and Landau theory: Consider a 
model of hard-core bosons with nearest-neighbor repul- 
sion on the triangular lattice, 



H=-J x £~ +h.c.)+ J*5> 
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(i) 



For large Jj_ / J z = A the ground state of this Hamilto- 
nian is a uniform superfluid. In the interaction dom- 
inated regime, for A < 0.2, quantum Monte Carlo 
simulations || have shown that the uniform superfluid 
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is unstable to density-wave modulations at wavevectors 
±Q = (±47r/3,0), leading to a supersolid. Through the 
standard mapping (S^ — b\. = b i: Sf — n~ 1/2), this 
Hamiltonian is equivalent to an S = 1/2 XXZ spin model 
with ferromagnetic in-plane interaction Jj_ and an anti- 
ferromagnetic out-of-plane interaction J z . Spin- wave cal- 
culations on this XXZ model |9], valid at large- £>, show 
that the phase transition into the supersolid is caused 
by the vanishing of the roton gap at ±Q, which leads 
to roton condensation. They also capture the precise 
structure of the supersolid phase of model (1) which is 
now known from direct numerics for 5 = 1/2 8]. Fur- 
ther, both techniques agree that the transition from the 
superfluid to the supersolid is continuous which is also 
indicated from a GL analysis outlined in Melko et al. Q . 
Since we are interested in the consequences of a contin- 
uous superfluid to supersolid transition, we focus on this 
specific model. We will use (semi) classical analyses of 
the XXZ model, supplemented by GL theory, as a reli- 
able guide to the physics of model (1). 

Keeping the relevant wavevectors q = 0,±Q to de- 
scribe the superfluid to supersolid transition, the super- 
fluid order, ip(r), and the deviation of the density from 
half-filling, Sn(r) = n(r) — 1/2, can be expressed as 
tp(r) ~ M r ) + V' Q (r)e iQ - r + ip_ Q (r) e - 4Q - r and Sn(r) ~ 
m(r) -|- p(r)e lC t' r + p*(r)e~ lC *' r . The complex numbers 
ipo,ip± Q ,P and the real number m are, thus, order pa- 
rameters in the GL theory 0, IToL Hlj . For the model 
(1), ipo is nonzero in, both, the superfluid and the super- 
solid, while ip± Q , m, p are nonzero only in the supersolid. 
Since we are interested in a hard-spin formulation here, 
it is convenient to work with the phase and amplitude 
of the superfluid order parameters ipo,ip ±Q . Terms such 
as ip^tp*_ Q ipo and V , qV'* q ' ( /'o m tnc GL functional lock 
the phases of these different components, so that ipo = 
|^o|e l¥> , ip ±Q = IV'iea l 6 *^ are determined by three ampli- 
tudes and a single phase variable <p. The amplitudes are 
constrained by the |V'o| 2 + |^ Q | 2 + |V- , _ Q | 2 + ™ 2 + 2|/9| 2 = 1, 
which fixes the spin length (to unity). 

The GL functional / = f p + f m + + f c , where 

f P = ap|p| 2 +ffp|Vp| 2 + Up |p| 4 + u; p Re(p 6 )... (2) 

fm = a m TO 2 +g m (Vm) 2 + u m m 4 . . . (3) 

U = Q S \V^\ 2 + ... (4) 

f c = -Aie 2 |p| 2 |V^| 2 + A 2 mRe(p 3 ) + A 3 [p 2 (V^ * + 

^qV'o) + cc] + A 4 [topV_ q V'o + c - c -] ■ ■ ■ ( 5 ) 

where we have only displayed terms relevant to our anal- 
ysis below. These terms respect the particle-hole and 
lattice symmetries of the Hamiltonian in Eq. (1), and 
could be, in principle, derived using functional integral 
methods as outlined for a related model in Ref. ^3] • Here 
we treat the coefficients as phenomenological parameters. 
In Eq. (J2J, ot p = a(A — A c0 ) measures the distance 
to the superfluid-supersolid transition point with a > 0. 
Spin-wave theory and Monte Carlo calculations indicate 



w p < 0, which prefers p to be real in the supersolid. In 
f c , the term Ai couples the modulated solid order and the 
superflow current QsS/ip. The A2,3,4 terms induce effective 
"magnetic fields" upon the different order parameters. 

Current driven collapse of the roton gap: For su- 

perfluids close to a supersolid instability, external per- 
turbations that suppresses the kinetic energy relative to 
the interaction energy can drive the system into the su- 
persolid phase. Consider a uniform superflow, say in 
the x direction, which introduces a phase difference be- 
tween neighboring sites, 6 — <j)j — 4>i- The phase gra- 
dient effectively suppresses the nearest-neighbor kinetic 
energy from A to Acos(0j — <j)j). Carrying through 
the calculation of spin-wave fluctuations around such 
a state with superflow, we find a spin-wave dispersion 

with a roton gap u>(Q) = l-5y 6A(A — 1/2), where 

A = A[2cos((5/2) + cos<5]/3, reduced by the current 
flow. The point where the roton gap collapses is thus 
shifted upward to A c ~ A c o[l + (<5/2) 2 ], where in the last 
step we assumed 5 <C 1. Recalling the current density 
J ~ sin S w S, we conclude the boost of A c is quadratic 
in the supercurrent density close to the transition, and 
the superflow induces supersolid ordering. 

In the GL approach, superflow-induced supersolid or- 
dering can be described phenomcnologically by the cou- 
pling term Ai between the supercurrent and the super- 
solid order parameter p in Eq.©. This term effectively 
shifts a p to a' p — a p — AiJ 2 , where the supercurrent 
J s = g s \Vtp\. We thus see that one can induce a transi- 
tion from a superfluid {a' > 0) into a supersolid (a' < 0) 
by a superflow which leads to a sufficiently large |Vy|. 
The "critical current density" required to obtain the su- 
persolid is J c = yja.p/\i. This is analogous to superflow 
induced spin-density modulation and gradient cou- 
plings inducing competing phases [14| considered earlier 
only within GL theory. We next turn to implications of 
this for a vortex in the superfluid phase. 

Vortex in the superfluid phase: For a superfluid 
vortex, the current density increases upon approaching 
the core as J s = Q s /r. The above result then indi- 
cates that the superflow would cause a supersolid pat- 
tern to be stabilized in a region with characteristic ra- 
dius R v — Q s ^J\\/a p around the core. It is clear from 
the GL functional that a p ~ E^ ot , since the roton gap de- 
termines the distance to the transition. We thus identify 
Rv ~ Qs/E m t. We have confirmed this by microscopic 
calculations where we have evaluated the roton gap, the 
superfluid stiffness D s (which has energy units and can 
be viewed as Q s /m* with an effective mass m*) and the 
"critical current" using spin-wave theory. 

To study, further, the interplay of the various order 
parameters near a vortex, we numerically studied model 
JU in the presence of an orbital magnetic field, by re- 
placing t in Eq. QJ with ty = texp(z j l A ■ dl). We 
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considered an L x L lattice with periodic boundary condi- 
tions, and a vector potential corresponding to a uniform 
magnetic field with the total flux chosen to be one flux 
quantum. The ground state then corresponds to a single 
vortex within the cell. Transforming to an effective spin 
model, and treating the spins as classical unit vectors, 
we found the ground state spin configuration {S r } us- 
ing a simulated annealing algorithm. The ground state 
corresponds to a single vortex with its core on a site of 
the lattice. The local order parameters at each r, i.e. 
p, m, ipo, and ip ±Q , were extracted from {S r } by coarse 
graining over six nearest neighbor sites. 

The inset to Fig. ^ shows the map of the boson den- 
sity Sn(r) around a vortex for A = 0.505 — it is clear 
from this that supersolid order is induced in the vicinity 
of the vortex core as expected from earlier arguments. 
Each data point in the plotted profiles in Fig. ^ repre- 
sents the angular average taken at a fixed distance from 
the core. It is apparent that as the superflow induces a 
nonzero p, the order parameter |-0o| is suppressed. At the 
same time, |"0±qI and to are generated but they are in- 
significant except very near the vortex core. We explain 
these features below using GL theory. 

Outside of a few lattice spacings from the vortex cen- 
ter m ~ 0, ^i ±p ~ and we will ignore them to be- 
gin with. Assuming p is small, and \ip\ is almost con- 
stant, the free energy density relevant to p reduces to 
/far « {a p -X ig 2 \V<p\ 2 )\p\ 2 +g p \Vp\ 2 +u p \p\ 4 + . . ., where 
we have dropped the higher order w p term but implicitly 
use w p < to consider only real p solutions below (the 
numerical solution for the vortex yields an almost purely 
real p). It is clear that for small r, \Vip\ ~ 1/r is large 
and renders the state with p — unstable. Minimiz- 
ing / with respect to p* , introducing s(r) — p(r)2u p /a p , 
one finds s(r) obeys the differential equation £z [d 2 s/dr 2 + 
(1/r) {ds/dr)] -s 3 + (R 2 Jr 2 -l) S = 0. Here ^ = Vfl^ 
is the coherence length for p, and R v — g s yj X\/ a p as 
defined before is the length scale associated with the su- 
persolid order. The solution to this equation at large r 
is s(r) ~ e~ r ^ p / \fr. This solution breaks down close 
to the vortex core as p saturates to a finite value at 
r = 0. The numerical result for p(r) can be fit well 
over the region r > 1 with p(r) — Ae^ 7 "^ 1 tanh(y^/r), 
as shown by the solid line in Fig. The length scale 
Rv ~ [( r ' 2 / 3 ( r ))/(p( r ))] 1 ^ 2 : where the (.) denotes the spa- 
tial average with the vortex at the origin. Once we have 
a nonzero \p\ 2 (far from the core) it suppresses, via the 
hard-spin constraint, the uniform superfluid order pa- 
rameter IV'ol 2 with [IV'ol 2 +2|p| 2 ] ~ 1 since m,ip ±Q ~ 0. 
Next we turn to these order parameters which are small 
everywhere. We see from Eq. 10 that X 2 Ke(p 3 ) acts as 
a magnetic field for to, while X^p 2 ^ induces ip Q . Thus 
we expect to ~ Re(p 3 ) and IV^qI ~ IpPIV'oI (except very 
close to the core), as borne out by our numerics. 

Tunneling: We next analyze electron tunneling into the 




FIG. 2: Tunneling DOS for an /-wave triplet SC on the tri- 
angular lattice. We set the spinon hopping tf = 1, the spinon 
pairing gap Af = 0.4, and set J± — J z ~ 0.08 in Eq. (1) 
for the condensate. Low energy roton modes do not influence 
the V-shaped nodal spectrum due to coherence factors. The 
lower energy peaks are elastic spinon peaks, secondary peaks 
arise from condensate fluctuations. At smaller S, the inelastic 
contribution is larger while the weaker condensate dispersion 
shifts the secondary peak to lower energy. 

uniform SC to see if it probes the low-energy roton ex- 
citation when the system is on the brink of becoming a 
supersolid. We use a slave boson formulation, writing 
the electron operator as c CT (r, r) = / CT (r, r)6(r, t) where 
fa- is a spinful fermion (spinon) and b is a charged boson. 
Assuming that gauge fluctuations and interactions are in- 
nocuous in the superconducting state, the electron Green 
function can be factorized as G%(r, t) = G£(r, r)G fc (r, t). 
In mean-field theory, when the boson is condensed, we 
can replace G b (r, t) ~ (b) 2 , so that G c ~ G' , their ratio 
being the condensate density (6) 2 . The electron tunnel- 
ing spectrum thus reflects the tunneling density of states 
(TDOS) associated with the spinon spectrum. 

Going beyond mean field theory, we assume that the 
bosons are described by an interacting Hamiltonian such 
as (I) and evaluate the boson Green function using spin- 
wave theory. Combining this with the spinon Green func- 
tion, for a model with spinon pairing, leads to the TDOS 

k 

+^ E ^ H 5 ^ + E * + Q «) + - E *~ "q)] ■ 

q#0,k 

Here u 2 ,, v 2 . are the usual superconducting coherence fac- 
tors associated with spinon pairing, while 7 2 is the co- 
herence factor associated with the Bogoliubov transfor- 
mation for diagonalizing the boson model 0| ■ The first 
term, which dominates at large-5, reflects the mean field 
result where the condensate does not fluctuate. The sec- 
ond term represents inelastic processes where tunneling 
inserts a spinon and excites condensate fluctuations. 

To show the effects of condensate fluctuations in a 
nodal SC on the triangular lattice, we plot in Fig. [21 



4 



the TDOS N c (uj) obtained for /-wave triplet pairing of 
spinons [16|. The two main features of the spectrum in 
Fig. J2J are: (i) At low energy the nodal spectrum is 
nearly unaffected by the condensate fluctuations, (ii) At 
higher energies, condensate fluctuations lead to a broad 
secondary peak. This arises from q regions where there 
is a large density of states for condensate excitations. 
These remain true even in the presence of low-energy ro- 
tons since the 7q coherence factor tends to reduce the 
contribution to the TDOS coming from near the roton 
wavevector Q. The energy difference between "elastic" 
and "inelastic" peaks scales with the superfluid stiffness. 

Experimental implications: Our results for the super- 
flow induced supersolid, the vortex core size scaling and 
the existence of secondary tunneling peaks rely mainly 
on the existence of a continuous superfluid-supersolid 
transition and a small superfluid density, largely inde- 
pendent of the detailed microscopies. Here we consider 
some implications for the cuprate superconductors. Tun- 
neling measurements on near-optimal superconducting 
Bi 2 Sr 2 CaCu 2 8+( 5 (T c = 89K) in a B = 5T magnetic 
field see a modulated TDOS around the vortex core [l7j . 
Although different explanations have been proposed for 
this pattern (e.g., Refs.JiO, El E3|)> it is possible that 
over the large region (with radius ~ 75A) in which the 
modulated pattern appears the system may be better 
viewed as a supersolid (which appears to be seen at low 
doping and B — in a related cuprate 0). If this super- 
solid has a charge modulation, the roton in this system 
at B — could be observed using inelastic X-ray scat- 
tering. We expect that the roton energy should scale as 
E ro t ~ D s /R v where D s is the bulk superfluid stiffness 
and R v is the radius of the region around the vortex core 
exhibiting a checkerboard modulation. Further, the ob- 
served modulation region around the vortex core should 
grow with underdoping upon approaching the supersolid. 
Both predictions could be tested experimentally. We 
have also shown that condensate fluctuations can lead to 
inelastic secondary peaks at higher energies, while leav- 
ing the low energy V-shaped spectrum unaffected. This 
is broadly consistent with tunneling data in the cuprates 
where regions with a small superfluid density and broad 
secondary tunneling peaks (so-called "large gap" regions) 
coexist with "small gap" regions of larger superfluid den- 
sity and single peaks |lj|. However, the fact that the high 
energy structure carries information about the ordering 
wavevector in the cuprates ^| does not appear to be 
a feature of the model we have studied. This may need 
coupling to disorder which can simultaneously lead to low 
superfluid density and supersolid order in certain regions 
and also to scattering of quasiparticles off the disorder 
induced supersolid. Further, one needs to address pos- 
sible bond-order instabilities [20I Eij . |24| . We hope our 
work will also stimulate experiments to study supersolids 
and vortices [2j| using cold atoms on optical lattices, and 



vortices in superfluid 4 He near freezing pressure. 
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